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Abstract: GRAPE-Dilepton is a Monte Carlo event generator for dilep- 
ton production in ep collisions. The cross-section calculation is based on 
the exact matrix elements in the electroweak theory at tree level. The 
dilepton productions via 77, 7Z , Z°Z° collisions and via photon internal 
conversion are taken into account. In addition, the effects of the Z° on/off- 
shell production are also included. The relevant Feynman amplitudes are 
generated by the automatic calculation system GRACE. The calculation of 
the proton vertex covers the whole kinematical region. This generator 
has an interface to PYTHIA and SOPHIA to obtain complete hadronic final 
states. 



PROGRAM SUMMARY 



Title of program: GRAPE-Dilepton (vl.l) 

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. 
Ireland and from http://www-zeus.desy.de/~abe/grape/. 

Operating system under which the program has been tested: UNIX 

Programming language used: Fortran77 

Memory required to execute with typical data: 7 Mwords for integrations, 9 Mwords 
for event generations 

Keywords: dilepton, lepton-pair, ep collision, Bethe-Heitler, Z boson, dipole form 
factor, hadron tensor, lepton tensor, structure function, parton density, GRACE 

Nature of physical problem: A precise estimation of the cross section of the elec- 
troweak dilepton production in ep collisions is required in various physics analyses, 
where 8~48 Feynman diagrams can contribute. 

Method of solution: The automatic calculation system GRACE is used to obtain all of 
the relevant helicity amplitudes. The phase space is divided into the 3 regions accord- 
ing to the kinematics at the proton vertex, and the 3 different calculation methods 
are applied. The radiative corrections are included using the structure function and 
the parton shower methods. 

Restrictions on the complexity of the problem: Higgs, the proton-Z° coupling and 
lepton pair production through photon radiation from the proton are not included. 
The contribution from the resolved photon, i. e. Drell-Yan process in ep collisions is 
not included. 

Typical running time: 1 hour for a cross-section integration and 1 msec per 1 event 
for an event generation 



1 Introduction 



In the study of electron/positron-proton (ep) collisions, a precise estimation of the 
dileptonP] production cross-sections in the electroweak (EW) interaction is important 
since it could become a significant background for various physics analyses such 
as, for example, exclusive J/t/i or T production and new physics searches. So far 
only the generator LPAIR |3J has been used in experimental analyses to estimate 
the dilepton background ||. The calculation of LPAIR is based on the diagrams 
of the photon-photon collision process ||, so-called two-photon Bethe-Heitler (2- 
7 BH), corresponding to the diagrams of Fig. (IJ-(a) or Fig.Q-(a) with the photon 
contribution only. This process is dominant in most of the phase space. It is, 
however, expected that in the region of low invariant masses of the dilepton system, 
QED-Compton type (CO) diagrams (i.e. photon internal conversion process) as seen 
in Fig.[I|-(b) and Fig. 0-(b) become dominant. In the high mass region, there is 
an additional interesting process, i.e. Z° production, which is implemented in the 
MC event generator EPVEC @|. EPVEC, however, does not include 2-7 BH nor CO 
diagrams. In the di-e channel, interference effects of the final state e~e~ or e + e + 
should also be taken into account, which are included neither in LPAIR nor in EPVEC. 

In this paper, a new MC event generator GRAPE-Dilepton for the dilepton pro- 
duction in ep collisions is presented. The FORTRAN code to calculate the Feynman 
amplitudes is generated by GRACE || which is an automatic calculation system. GRACE 
has been used mainly for e + e~ interactions so far. This is the first time for GRACE to 
be applied to the case where there is a composite particle (i.e. proton) in the initial 
state. GRAPE stands for a GRA ce-based generator for Proton- Electron collisions. 

This generator has the following features. 

• The cross-section calculation is based on the exact matrix elements in the elec- 
troweak theory at tree level. Not only 2-7 BH but also the dilepton productions 
via 7Z and Z°Z° collisions are taken into account. CO and Z° on/off-shell pro- 
duction are also included. Interference effects of the final state e ± e ± are taken 
into account in the di-e channel. It is possible to select any sub-set of diagrams 
in the calculation. 

• All fermion masses are kept non-zero both in the matrix elements and in the 
kinematics, which makes it possible to use this program with arbitrary small 

lr The word dilepton represents di-electron(di-e), di-muon(di-/i) and di-tau(di-r) in this paper. 
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scattering angles of and/or small invariant masses of dilepton down to the 
kinematical limits. 

• The calculation of the proton vertex covers the whole kinematical region by 
dividing it into 3 categories of elastic, quasi-elastic and DIS (Deep Inelastic eq 
Scattering) processes, as described in the next section in detail. 

• Both of Initial State Radiation (ISR) and Final State Radiation (FSR) can be 
included. 

2 Physics aspects 

This generator simulates the ep interaction: e^ n ^ p {m) — > e ± l + l~X where and 
p {in) indicate the electron/positron and the proton in the initial state respectively, 
and l + l~ are the scattered electron/positron and the produced dilepton respectively. 
The relevant processes are classified into 3 categories using the negative momentum 
transfer squared at the proton vertex (Qp) and the invariant mass of the hadronic 
system (M had ); 

Ql = -{p^-(p e ±+p l+ +p l -)} 2 > (!) 
K* = {b e±(in) +p p ^) - (p e± +p l+ +?v)} 2 > ( 2 ) 

where p , and p ,. , are the 4- momenta of the incoming lepton and the proton after 
ISR, respectively. p e± and p i± are those of the scattered lepton and the produced 
leptons before FSR, respectively. The 3 categories are 

• Mhad = M p (elastic), 

• Q 2 P < Qmin 0R M p + M tt° < M had < M cut (quasi- elastic) , 

• Ql > Qmin AND M had > M cut (DIS), 

where M p and M^o are the masses of the proton and the neutral pion, respectively. 
Qmin is set to around 1 GeV depending on the Parton Density Function (PDF) used 
in the DIS process. The recommended value for M cut is 5 GeV. 

For the elastic process, the diagrams in Fig.[I| are calculated with the following 
dipole form factor for the proton-proton-photon vertex (r^ p7 ) with the on-shell pro- 
ton. The general form of the elastic proton vertex can be written as 

T»p= e p (f^QD^ + l^fFiiQl) W<b) (3) 
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where e p indicates the electric charge of the proton, q is the 4-momentum transfer 
at the proton vertex (q 2 = —Qp), Fi(Q 2 ) and F 2 (Q p ) are the 2 independent form 
factors, and k p is the anomalous magnetic moment of the proton (see, for example, 
||.). The electric and magnetic form factors G P E {Q 2 ) and G P M {Q 2 ), respectively are 
defined as follows, 




F 1 (QD + k p F 2 {QD 



(4) 



Using the Gordon decomposition and the scaling law of the form factor, 

G P E (Ql)=G p M (Q 2 p )/M ( 5 ) 
the following formula which is used in this program is obtained, 

V P 1 + 4A^ / 

where fi p = (1 + K p )[i B) fi B is the Bohr magneton, and p p ( ou t) indicates the 4- 
momentum of the scattered proton. G E (Q p ) is calculated according to the formula 
of the dipole fit, 



0.71 GeV\ 

The only difference between the elastic and the quasi-elastic processes is the treat- 
ment of the proton vertex and the simulation of the hadronic final state. The quasi- 
elastic proton vertex can be described using the hadron tensor in the following form 
assuming parity and current conservation (for example, see 0.), 

= W x (-g^ + ^X- 
\ T 

_i_ W \ ( _ ^W") 1 n n } ( V _ Pp(in) Q i 

2 Ml V Pp(in) ~l^ q ) \ Ppiin) ~l^ q 

W±(Qp, Mhad) and W2(Q p , Mhad) are the electromagnetic proton structure functions. 
The hadron tensor is contracted with the lepton tensor L^ u numerically to obtain 
the cross section, 

da ~ L,„r. (9) 
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In this version, W\ and Wi are parameterized with Brasse et al. |7j for Mhad < 2 GeV 
(the proton resonance region), and with ALLM97 || for Mhad > 2 GeV. These two 
parameterizations are based on fits to the experimental data on the measurement of 
the total 7*]? cross-sections. The exclusive hadronic final state is generated using the 
MC event generator SOPHIA || in the event generation step. 

For the DIS process with the Quark Parton Model, the diagrams in Fig.|2| are 



calculated. PDFLIB |10[ is linked to obtain parton densities with Q p as a QCD 
scale. The simulation of the proton remnant and the hadronization are performed 
by PYTHIA ||11J| . It should be noted that the lowest order calculation in this process 



is valid only for the region of 

u= l^-^+^J} 2 ! > 25GeV 2 , (10) 

where p q(in) is the 4-momentum of the incoming quark. The value of u corresponds 
to the virtuality of the u-channel quark in the diagrams in Fig.Q-(b),(c). When it 
is nearly or smaller than 25 GeV 2 , the lowest order calculation is not correct as 
explained in 0] since QCD corrections become large. In this case, the dilepton pro- 
duction should be treated as Drell-Yan process between the proton and the resolved 
photon from the beam lepton, which is not implemented in this program. The cut: 
u > 25 GeV 2 is explicitly applied in this program if the diagrams other than BH are 
included. 

The effect of ISR is included in the cross-section calculation using the structure 
function method described in fll2"f, where the momentum transfer squared on the 



beam lepton, i.e. \p , — p , } 2 is used as a QED scale. When ISR turns on, 
the correction for the photon self energy, i.e. the vacuum polarization, is included 
according to the parameterization in [H| by modifying photon propagators. FSR is 



performed by PYTHIA using the parton shower method when the event is generated. 
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(a) Bethe-Heitler type diagrams 




(c) Z on/off-shell production 



Fig. 1: Feynman diagrams included in the (quasi-)elastic process. e={e + ,e }, 
\ ± —{e ± ,/i ± ,T ± }. N means a (dissociated) proton or a nucleon resonance. 
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(b) QED-Compton type diagrams 



Fig. 2: Feynman diagrams included in the DIS process. e={e + , e }. l ± ={e ± , t^} 
and q={tfff<i§' > <fiW}. 
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3 Program structure 



Physics events are generated with the 2 steps; the MC integration step by the exe- 
cutable: integ and the event generation step by the executable: spring, as illustrated 
in Fig.|3|. In both steps, the program is controlled by an ASCII file: grape . cards. 
The file is read by the executables with help of FFREAD [14]. The contents of 
grape . cards are explained in the next section. 

In the integration step by the executable: integ, an effective total cross-section 
(in unit of pb) and probability distributions are calculated by BASES fl5fl . The results 
are stored in a file: bases .rz which has the Ntuple format provided by the HBOOK 
package [IB]. At the same time, the information related to the convergency status 
of the integration is output into an ASCII file: bases .result. 

In the event generation step by the executable: spring, unweighted events are 
generated. This is done by an routine: SPRING [13] according to the probability dis- 
tributions in bases .rz. The results of the event generation are stored in the PYTHIA 
common block /PYJETS/. After filling /PYJETS/, spring calls a routine: USRSTR in 
which user specific procedures are put. Its template is found in the appendix. The 
event information in /PYJETS/ is also available in a Ntuple file: grp. rz. 

The calculated cross-section is found in bases .result or at the end of the stan- 
dard output from spring. The status of the event generation is output into an 
ASCII file: spring. result. Looking at the file, users should find a reasonable agree- 
ment between generated distributions by spring and calculated ones by integ. The 
procedure to make the executables is described in the README file. 



4 Input data cards 

The input data in grape . cards are explained in this section. All of the items 
are optional and are set to default values if not specified. Default values are written 
in the brackets starting with D=. The items are not explicitly displayed in case that 
they are the only one for their cards. 

• KFLBEAM 

KF code of the lepton beam (INTEGER, D=-ll); ll:electron, -ll:positron. 

• EPOL P 9 cj) 

Polarization of the lepton beam (REAL); 
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grape.cards 



1st step 



|integ| 



1/ 



bases, rz 



bases. result 



2nd step — > ; - |spring| - 



Generated events 
Elastic Quasi-elastic DIS 

^ r 



(SOPHIA) (PYTHIA) 



(USRSTR)^- /PYJETS/ 



1/ 



V 



v;.. 



spring.result 



grp.rz 



I I Executable 



Input/Output file 



*.rz Ntuple file 



*. result ASCII file 



Fig. 3: Flowchart for the program structure 



P = degree of the polarization in the range [—1, 1] (D=0.), 
9 = polar angle of the polarization vector in degree (D=0.), 
= azimuthal angle of the polarization vector in degree (D=0 . ) . 

The positive direction of the z-axis on the polarization vector is in the direction of 

the lepton beam. 

• EBEAM 

Lepton beam momentum in MeV/c (REAL, D=27520.). 

• PBEAM 

Proton beam momentum in MeV/c (REAL, D=820000.). 

• PROCESS 

Process type of the proton vertex (INTEGER, D=l); 
lrelastic, 2:quasi-elastic, 3:DIS. 

• LPAIR 

Dilepton channel (INTEGER, D=2); l:di-e, 2:di-//, 3:di-r. 

• ISR 

Initial state radiation flag for the beam lepton (INTEGER, D=l); 0:off, l:on. 



8 





IVlJi/Hljrrj 


(Quarks 




LlhjU/hj W/CU/Z 


1 


1234 


u + u + d + d 


Yes 


No 


1 


123456 


u+u+d+d+s+s 


Yes 


No 


1 


12345678 


u+u+d+d+s+s+c+c 


Yes 


No 


1 


1234567890 


u+u+d+d+s+s+c+c+b+b 


Yes 


No 


1 


17 


u + c 


Yes 


Yes 


2 


28 


u + c 


Yes 


Yes 


3 


35 


d+ s 


Yes 


Yes 


4 


46 


d + s 


Yes 


Yes 


3 


359 


d + s + b 


Yes 


Yes 


4 


460 


d + s + b 


Yes 


Yes 



Table 1: Possible combinations of QFLV and MERGE 

QFLV 

Scattered quark in the DIS process (INTEGER, D=l); 

l:u, 2:u, 3:d, 4:d, 5:s, 6:s, 7:c, 8:c, 9:6, 10:6, ll:t, 12:1 
MERGE 

Merging mode in the DIS process (INTEGER, D=0); 0:off. 

In some cases, contributions from different quarks can be included in the cross- 
section calculation adding the parton densities if the mass difference is negligible. 
The possible combinations of QFLV and MERGE are written in Table|l|. The 
mass of the quark specified with QFLV is used in the amplitude and the kinematics 
calculations. 
NGROUP 

Author group described in the PDFLIB manual (INTEGER, D=5). 
NSET 

PDF set described in the PDFLIB manual (INTEGER, D=5). 

The default is GRV94(LO). 

GRASEL 

Feynman diagram selection (INTEGER, D=3); 

1 : 2-7 Bethe- Heftier (without e ± e ± interference in case of di-e), 

2 : 2-7 Bethe- Heitler (including e ± e ± interference in case of di-e), 

3 : QED diagrams (i.e. all the diagrams except for the Z° contribution), 

4 : EW diagrams (i.e. all the diagrams), 
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13 : QED-Compton type diagrams only, 

14 : Z° production diagrams only. 
In case of di-/x, r, the first and the second selections give the same result. 
ITMX1 

Number of iterations in the grid optimization step of BASES. (INTEGER, D=4). 

This should be larger than 2. 

ITMX2 

Number of iterations in the integration step of BASES (INTEGER, D=10). 

This should be larger than 5. 

NCALL 

Number of sampling points in each iteration of BASES (INTEGER, D=1000000). 
This should be large so that any accuracy of each iteration in the integration step 
of BASES is better than 0.5 %. 
NGEN 

Number of events to be generated by spring (INTEGER, D=100). 
NMOD N mod 

Printing a message per N mo d events in the event generation (INTEGER, D=1000). 
PSISR 

Switch for the initial state parton shower by PYTHIA (INTEGER, D=l); 0:off, Iron. 
This has an effect only on event generations of the DIS process. No effect on elastic 
and quasi-elastic events. This item is copied to MSTP(61) in the PYTHIA common 
block /PYPARS/. 
PSFSR 

Switch for the final state parton shower by PYTHIA (INTEGER, D=l); 0:off, l:on. 
This item is copied to MSTP(71) in the PYTHIA common block /PYPARS/. 
PSBRA 

Parton shower branchings in PYTHIA (INTEGER, D=2); 

1:QCD, 2:QCD+QED. 
This item is copied to MSTJ(41) in the PYTHIA common block /PYDAT1/. 
PSSUP 

Suppression of the PYTHIA parton shower (INTEGER, D=0); 0:off, >=l:on. 
This item is copied to MSTJ(40) in the PYTHIA common block /PYDAT1/. 
PYDECAY 

Switch for fragmentation and decay in PYTHIA (INTEGER, D=l); 0:off, Iron. 
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No effect on efastic and quasi-elastic events. 

This item is copied to MSTP(lll) in the PYTHIA common block /PYPARS/. 
PRIPT 

Primordial k t distribution in the proton (INTEGER, D=l); 

0:off, l:gaussian, 2:exponential. 
No effect on elastic and quasi-elastic events. This item is copied to MSTP(91) in 
the PYTHIA common block /PYPARS/. 
PYLIST 

Printing the contents of /PYJETS/ (LOGICAL, D=TRUE). 
NLIST 

Number of events whose /PYJETS/ is printed out (INTEGER, D=10). 
NTPYT 

Output of generated events into a Ntuple file: grp. rz 
from the PYTHIA common block /PYJETS/ (LOGICAL, D=FALSE). 
The meanings of the Ntuple variables are in the following, 
npy : Number of particles (integer) 

px(l:npy) , py(l:npy) , pz(l:npy) : x,y,z-component of momentum 

in GeV/c (real*4) 

pe(l:npy) : Energy in GeV (real*4) 
pm(l:npy) : Mass in GeV (real*4) 
kf (l:npy) : KF code (integer) 
sta(l:npy) : Status code (integer) 

mot (1: npy) : Line number of the mother particle (integer) 
Q2RNGME Mm Max 

Range for the negative momentum transfer squared at the electron vertex Q\ 
without ISR (REAL), i.e. Q 2 e = —{p e±(in) — P e ±} 2 where P e±(in) is a 4-momentum of 
the incoming lepton after ISR. 

Min = the minimum in GeV 2 (D=0.), 

Max = the maximum in GeV 2 (D=1.E20). 
In case of di-e with e ± e ± interference, smaller one of the two Q 2 e values is used. 
Q2RNGOB Min Max 

Range for the negative momentum transfer squared at the electron vertex Q\ 
including ISR (REAL), i.e. Q\ = — {P e±{m) ~ P e ±} 2 where P e±{in) is a 4-momentum 
of the incoming lepton before ISR. 
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Min = the minimum in GeV 2 (D=0.), 

Max = the maximum in GeV 2 (D=1.E20). 
In case of di-e with e ± e ± interference, smaller one of the two Q\ values is used. 
MHAD Min Max 

Range for the mass of the hadronic system M had (REAL); 

Min = the minimum in GeV (D=1.08), 

Max = the maximum in GeV (D=1.E20). 
No effect on elastic events. 
Q2P Min Max 

Range for the negative momentum transfer squared at the proton vertex Q 2 (REAL) ; 

Min = the minimum in GeV 2 (D=0 . ) , 

Max = the maximum in GeV 2 (D=1.E20). 
In case of the DIS process, Q 2 is used as a QCD scale for PDF. 



THMIN 


0(1) 

min 


0(2) 
min 


0(3) 0(4) 
min min 


THMAX 


0(1) 
"max 


fl(2) 
Umax 


fl(3) fl (4) 
Umax Umax 


EMIN 


E {1) 

mm 


E {2) 

mm 


E i3) E {A) 

mm mm 


EMAX 


J-Jmax 


P (2) 
J-*max 


E {3) E {4) 

J-'max -'-'max 


PMIN 


p(l) 

min 


p(2) 
min 


p(3) p(4) 

min min 


PMAX 


p(l) 
±max 


p(2) 
max 


p(3) p(4) 
max J^max 


PTMIN 


Pt W 

min 


Pf m P/ (3) p (4) 

'"min '"min '"min 


PTMAX 


Pt {1) 


p f (2) p (3) p ,(4) 
J Lmax -L Lmax ± ^max 



for scattered proton or quark, 
for scattered e ± , 
for produced l T , 
for produced / ± . 

The above 8 data cards are used for describing the detector cut in the laboratory 
frame (REAL) . Each final state particle is required to satisfy the following, 



(1) 
(2) 
(3) 
(4) 



u min ° ^ °max 



l(0 
m 

>(0 



AND E m \ n <E<El\ x 



7(0 
nir 

t (0 



f (0 



AND P^ n <P<P^' a x AND Pt K Z n < Pt < Pt'ma, 

where ^(degree), E'(GeV), P(GeV/c) and Pt(GeV/c) indicate polar angle, energy, 
momentum and transverse momentum, respectively. The default values correspond 
to not applying this cut. 
THPTMCT min 

max 
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PTMXCT Pt min 

Using the above 2 data cards, final state leptons are required to satisfy the follow- 
ing (REAL), 

@min <^ Q < 9'rnax AND Pt m i n < Pt < Pt max 

where Pt M (GeV/c) indicates the maximum transverse momentum among the 3 

final state leptons (e ± , l T , l^), and ^(degree) is the polar angle of the lepton with 

Pt M . The default values correspond to not applying this cut. 

MASSLL Mini Maxl Min2 Max2 

Range for the mass of the produced dilepton system (REAL); 

Mini = the minimum in GeV (D=0.), 

Maxl = the maximum in GeV (D=1.E20). 
In case of di-/i, r, Mini and Maxl are not used. 
In case of di-e with e ± e ± interference, there are two masses; 

' ' e^e ' e'e~ 

(Mg+ < Mg+ ), and they are required to satisfy the following, 

Mini < m£ ] < Maxl AND Min2 < M^ ] < Max2. 
MASSELL Min Max 

Range for the mass of the final state lepton system of e ± l T (REAL); 

Min = the minimum in GeV (D=l.), 

Max = the maximum in GeV (D=1.E20). 
MASSQLL Min Max 

Range for the mass of the scattered quark and produced dilepton system of q l + l~ 
(REAL); 

Min = the minimum in GeV (D=5.), 
Max = the maximum in GeV (D=1.E20). 
This cut has an effect only on the DIS process. 

In case of di-e with e ± e ± interference, smaller one of the 2 values is used. 
IVISI N msi 

thevmin ei\ n etl etl e { ± 

THEVMAX 9 m L 9 max dmax Omax 
EVMIN E {1) E {2) E i3) E W 



mm mm mm mm 



PVMAX F (2) F (3) F (4) 

■ j v ivx^-v^v J^max J^max J^max J-^max 



PTVMIN PC n Pt% m P^L Pt { ± 

PTVMAY P-f^ PV-( 4 ) 

■ ■ v iv±^-vyv l Umax J- l^max J- l^max J- l<ma 

(1) : for scattered proton or quark, 
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(2) : for scattered e ± , 

(3) : for produced l T , 

(4) : for produced l ± . 

The above 6 data cards are used for describing the detector cut in the laboratory 
frame (REAL except for N visi : INTEGER). N visi particle(s) are(is) required to satisfy 
the following, 

€L < e < 9 { l x AND E± <E< E$ ax AND Pt% in < Pt < Pt% ax 
where ^(degree), E(GeV) and Pt(GeV/c) indicate polar angle, energy and trans- 
verse momentum, respectively. As for N visi , D=-l, which corresponds to not apply- 
ing this cut. The test run at the end of this paper is instructive for understanding 
this cut. 



5 Summary 

A new Monte Carlo generator for dilepton production in the framework of the 
electroweak theory was presented. The whole kinematical region on the proton vertex 
is covered. This generator can be used for quantitative and precise estimations of 
processes which come in addition to the two-photon Bethe-Heitler contributions. 
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Appendix A. Event store in /PYJETS/ 

Each line in the PYTHIA event store has the following meaning according to the 
PYTHIA convention. 



Line number 


Meaning 


1,2 


Beam particles (l:p, 2:e ± ) 


3,4 


Partons from the beam particles before ISR 
(3:from p, 4:from e^) 


5,6 


Partons from the beam particles after ISR; 

the initial state in the matrix element calculation 
(5:from p, 6:from e^) 


7,8,9,10 


Final state particles before FSR; 

the final state in the matrix element calculation 
(7:p, 8:e ± , 9:P, 10:/*) 


11~ 


Final state particles after FSR, fragmentation and decay 



In case of the (quasi-)elastic process, 

• a parton from the beam proton is the proton itself, i.e. the 1st and the 3rd 
lines are the same, 

• the beam proton always makes no ISR, so that the 1st, 3rd and 5th lines are 
the same, 

• the 2nd and the 4th lines are also the same. 

In case of the DIS process with both ISR and FSR, all of the lines have different 
contents in general. 

As for di-e events, there are 2 identical particles in the final state. In GRAPE, those 
2 particles are distinguished in the following way; 

• in case of e ± e ± interference off, a particle stored in the 8th line is a scattered 
lepton, and one in the 10th line comes from the 2-7 collision (i.e. a produced 
lepton), 

• in case of e ± e ± interference on, a lepton stored in the 8th line has smaller 
transverse momentum than that of a lepton in the 10th line. 
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Appendix B. Routines/function related to users 



• Function DRN(ISEED) 

provides uniform random numbers. All other random number routines are linked to 
this one. This routine is stored in the BASES library (libbases . a). 

• Subroutine DRNSET(ISEED) 

performs an initialization for DRN(ISEED). This is also stored in libbases. a. 

• Subroutine READ_CARDS(LUN, filename) 
reads input data cards from grape . cards. 

• Subroutine SETMAS 

provides masses/widths of particles and the QED coupling constant. 

• Subroutine USRSTR(Ievt ,Ngen) 

can be modified by users to access the information on generated events. User initial- 
ization/termination procedures are also put in this routine. 
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Appendix C. User event storing routine: USRSTR 



This routine is called (N gen + 2) times where N gen is the number of generated events. 
The 1st call is for the user initialization and the last one for the user termination 
phase. The following is the template file prepared in the GRAPE package (usrstr.f). 

subroutine USRSTR (Ievt , Ngen) 
implicit NONE 



* Arguments 

integer Ievt, Ngen 

* Ngen : # of events to be generated 

* Ievt : Counter < 1 ===> Initialization phase 

* 1 - Ngen ===> Event generation phase 

* > Ngen ===> Termination phase 
* 

* PYTHIA common 



integer N , NPAD , K (4000, 5) 

double precision P(4000,5), V(4000,5) 
common /PYJETS/ N,NPAD,K,P,V ! ! ! Event Record ! ! ! 

integer MINT (400) 

double precision VINT (400) 
common /PYINT1/ MINT, VINT 

* (See PYTHIA manual for details.) 

* 

* Local variables 

integer LUN1, LUN2, LUN3 
parameter (LUN1=41, LUN2=42, LUN3=43) 

* (You can use the above logical unit numbers.) 
* 

******** Initialization of USER Event Storing ******* 
if (Ievt .LT. 1) then 

endif 

************* <<< USER Event Storing >>> ************ 
if ((Ievt .GE. 1). and. (Ievt . LE. Ngen)) then 

endif 

********* Termination of USER Event Storing ********* 
if (Ievt .GT. Ngen) then 

endif 

return 
end 
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TEST RUN INPUT AND OUTPUT 

An example of the input data cards 

One example is presented with the following condition; 

• process : e+ q -> e+ q {q={v)ld ] } with GRV94(LO) 0), 

• BH diagrams only, 

• 70 % polarization of the e + beam in the direction of the proton beam, 

• cuts: (1) & (2) & (3) & (4), 

(1) Q2 >1GeV 2 & M had >5GeV, 

(2) for scattered q, 9 > 10° & P t > 15GeV/c, 

(3) (invariant mass of [i + > 4GeV, 

(4) 2 of the following 3 requirements for the final state leptons are satisfied, 

- for e + : 5° < 9 < 175° & P t > 5 GeV, 

- for fj+ : 20° < 9 < 160° & P t > 3 GeV, 

- for : 20° < 9 < 160° & P* > 3 GeV. 

Following is the corresponding grape . cards. The same file is put in the directory 
sample. 

LIST 

NCALL 1400000 

C ====================================================================== 

C « Polarization of the Lepton Beam » 
C (1) (2) (3) 

EP0L -0.7 0. 0. 

C ====================================================================== 

C « Process in the Proton Vertex >> (l:elastic, 2 : quasi-elastic , 3:DIS) 
PROCESS 3 

C ====================================================================== 

C « Produced Lepton-pair >> (l:di-e, 2:di-mu, 3:di-tau) 
LPAIR 2 

C ====================================================================== 

C « Scattered Quark in DIS » 

C (l:u, 2:u-bar, 3:d, 4:d-bar, 5:s, 6:s-bar, 7:c, ... , 12:t-bar) 
QFLV 1 
MERGE 123456 

C ====================================================================== 

C « PDF set in DIS » (See PDFLIB manual.) 
NGR0UP 5 
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NSET 5 

C ====================================================================== 

c ###################################################################### 

c ====================================================================== 

C << Electroweak Dilepton Production » 
GRASEL 2 

C ====================================================================== 

C ###################################################################### 

c ====================================================================== 

C << Mass Range for the Hadronic System >> (only for quasi-elastic and DIS) 
MHAD 5. 300. 

C 

C « Q2 Range for the Proton Vertex » 
Q2P ~ 1. 1.E20 
C 

C << Cuts for each Final-state Particle >> 
C ^/q? <e+-> <!-+> <!+-> 



THMIN 


10. 




0. 




0. 




0. 




THMAX 


180. 


180. 


180. 




180, 




EMIN 


0. 




0. 




0. 




0. 




EMAX 


1. 


,E20 


1. 


,E20 


1. 


.E20 


1. 


E20 


PMIN 


0. 




0. 




0. 




0. 




PMAX 


1. 


.E20 


1. 


,E20 


1, 


.E20 


1. 


,E20 


PTMIN 


15. 




0. 




0. 




0. 




PTMAX 


1. 


.E20 


1. 


.E20 


1. 


.E20 


1. 


,E20 



c 

C « Mass cuts » 
MASSLL 4. 1.E20 

0. 1.E20 

C 

C << Cuts for One or Some of the Final-state Particles >> 
IVISI 2 



c 


<p/q> 


<e+- 


-> <l-+> <1+- 


THEVMIN 


0. 


5. 


20. 20. 


THEVMAX 


0. 


175. 


160. 160. 


EVMIN 


0. 


0. 


0. 0. 


EVMAX 


0. 


1. 


,E20 1.E20 1.: 


PTVMIN 


0. 


5. 


3. 3. 


PTVMAX 


0. 


1. 


.E20 1.E20 1. 



c === 
STOP 
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A part of the standard output from integ 



#################################################################### 
##**++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++**## 



##**++ 
##**++ 
##**++ 
##**++ 
##**++ 
##**++ 
##**++ 
##**++ 
##**++ 
##**++ 
##**++ 
##**++ 
##**++ 
##**++ 
##**++ 



GGGGGG 



RRRRRR 



G 
G 
G 
G 
G 



G 



GGGGGG 
G G 
G G 
GGGGG G 



R 
R 
R 
R 
R 
R 



RRRR 
R 

R 



A 
A A 
A A 
A A 
AAAAAAA 



RR 



A 
A 



PPPPPPP 

P P 

P P 

PPPPPPPP 

P 

P 

P 



EEEEEE 



GRAce-based generator for Proton-Electron collisions 
GRAPE-Dilepton_versionl . 1 



++**## 
++**## 
++**## 
++**## 
++**## 
++**## 
++**## 
++**## 
++**## 
++**## 
++**## 
++**## 
++**## 
++**## 
++**## 



Mar. 14 2000 

Comments/bug-report to Tetsuo ABE(tabe@post .kek.jp) 

(abeOmail . desy . de) 
##**++++++++++++++++++++++++++++++++++++++++++++++++++++ 

#################################################################### 



««««« This is an INTEGRATION step. »»»»» 



> DIS process 

(Scattering of e and u quark) 
> Muon-pair production 

grace 2.1(5) 

(c) Copyright 1990-1998 Minami-Tateya Group (Japan) 
»> Graph selection 



jselg 




1 





2 





3 





4 





5 





6 





7 





8 





9 


1 


10 





11 





12 





13 


1 


14 





15 





16 





17 





18 





19 





20 





21 





22 





23 





24 





25 
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: > Start of Kinematics Initialization 



* * * * * * * * * * Information (in Lab. frame) ********** 

(in unit of GeV) 



P of electrons 




27.5200 


P of protons 




820.000 


Mass of electron 




5.10999E-04 


Mass of proton 




0.938272 


sqrt(S) 




300.444 


P of CMS 




792.480 


E of CMS 




847.521 


gamma of CMS 




2.82089 


beta*gamma of CMS 




2.63770 



« Mass range for the hadronic system » 
Min. = 5.00000 GeV 
Max. = 300.000 GeV 

> PDFLIB Initialization started 

***** PDFLIB Version: 7.09 Released on 970702 at 16.05 
in the CERN Computer Program Library W5051 ***** 
***** Library compiled on 970702 at 16.05 ***** 



Parm = Nptype Ngroup Nset 
Val = 1.0000 5.0000 5.0000 

Nptype = 1 Ngroup = 5 Nset = 5 Name = "GRV94-LD" CrMode = -1 

Nfl = -5, LD = 1, Tmas = 180.00 GeV/c**2 

QCDL4 = 0.2000 GeV, QCDL5 = 0.1530 GeV 

Xmin = 0.10E-05, Xmax = 0.99999E+00, 

Q2min = 0.400 (GeV/c)**2, Q2max = 0.10E+07 (GeV/c)**2 



> PDFLIB Initialization finished 

> ISR for incoming lepton using Structure Function method 

========> End of Kinematics Initialization 

>>> e+ beam 

Date: 0/ 4/12 21:29 



* 












* 


* 


BBBBBBB 


AAAA 


SSSSSS 




SSSSSS 


* 


* 


BB BB 


AA AA 


SS SS 


EE 


SS SS 


* 


* 


BB BB 


AA AA 


SS 


EE 


SS 


* 


* 


BBBBBBB 


AAAAAAAA 


SSSSSS 


EEEEEE 


SSSSSS 


* 


* 


BB BB 


AA AA 


ss 


EE 


SS 


* 


* 


BB BB 


AA AA 


ss ss 


EE 


SS ss 


* 


* 


BBBB BB 


AA AA 


SSSSSS 


EEEEEE 


SSSSSS 


* 


* 
* 




BASES 


Version 5. 


1 




* 
* 


* 


coded by S.Kawabata KEK, 


March 


1994 


* 



« Parameters for BASES » 



23 



(1) Dimensions of integration etc. 



# of dimensions 

# of Wilds 

# of sample points 

# of subregions 

# of regions 

# of Hypercubes 



Ndim = 9 ( 50 at max.) 

Nwild = 7 ( 15 at max.) 

Ncall = 1399680 (real) 1400000 (given) 

Ng = 48 / variable 

Nregion = 6 / variable 

Ncube = 279936 



(2) About the integration variables 





-+ 




-+- 






-+- 




-+ 




i 




XL(i) 






XU(i) 




IG(i) 




Wild 


1 


-+ 




000000E+00 


-+- 


1 


000000E+00 


•+- 


1 


-+ 


yes 


2 





000000E+00 




1 


000000E+00 




1 




yes 


3 





000000E+00 




1 


000000E+00 




1 




yes 


4 





000000E+00 




1 


000000E+00 




1 




yes 


5 





000000E+00 




1 


000000E+00 




1 




yes 


6 





000000E+00 




1 


000000E+00 




1 




yes 


7 





000000E+00 




1 


000000E+00 




1 




yes 


8 





000000E+00 




1 


000000E+00 









no 


9 





000000E+00 




1 


000000E+00 









no 




•+ 




-+- 






■+- 




-+ 





(3) Parameters for the grid optimization step 
Max.# of iterations: ITMX1 = 4 
Expected accuracy : Accl = 0.2000 7« 

(4) Parameters for the integration step 
Max.# of iterations: ITMX2 = 10 
Expected accuracy : Acc2 = 0.0100 °/ 

Date: 0/ 4/12 21:29 
Convergency Behavior for the Grid Optimization Step 



<- Result of each iteration -> <- Cumulative Result -> < CPU time > 
IT Eff R_Neg Estimate Acc 7. Estimate (+- Error ) order Acc 7, ( H: M: Sec ) 



1 2 0.00 4.586E-02 6.435 4.586406(+-0.295122)E-02 6.435 0: 1:37.57 

2 29 0.00 5.205E-02 0.295 5.203249(+-0.015359)E-02 0.295 0: 7:23.86 

3 37 0.00 5.196E-02 0.201 5. 198557(+-0.008625)E-02 0.166 0:13:52.62 



Date: 0/ 4/12 21:29 
Convergency Behavior for the Integration Step 



<- Result of each iteration -> <- Cumulative Result -> < CPU time > 

IT Eff R_Neg Estimate Acc 7, Estimate (+- Error ) order Acc 7. ( H: M: Sec ) 



1 


38 





00 


5 


196E- 


•02 





194 


5 


196162(+- 





010077)E-02 





194 





20 


27 


03 


2 


39 





00 


5 


201E- 


•02 





198 


5 


198416(+- 


-0 


007197)E-02 





138 





27 





89 


3 


38 





00 


5 


201E- 


02 





395 


5 


198705(+- 





006793)E-02 





131 





33 


34 


96 


4 


39 





00 


5 


199E- 


-02 





191 


5 


198759(+- 





005605)E-02 





108 





40 


9 


71 


5 


38 





00 


5 


200E- 


•02 





195 


5 


198989(+- 





004905)E-02 





094 





46 


43 


94 


6 


39 





00 


5 


200E- 


-02 





193 


5 


199218(+- 





004405) E-02 





085 





53 


18 


60 


7 


38 





00 


5 


187E- 


02 





193 


5 


197301 (+- 





004032)E-02 





078 





59 


53 


16 


8 


38 





00 


5 


185E- 


-02 





193 


5 


195602(+- 





003739)E-02 





072 


1 


6 


27 


22 


9 


38 





00 


5 


205E- 


02 





193 


5 


196769(+- 





003505)E-02 





067 


1 


13 


1 


21 


10 


38 





00 


5 


191E- 


-02 





196 


5 


196186(+- 





003313)E-02 





064 


1 


19 


34 


81 



****** end op BASES ********* 
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<< Computing Time Information >> 



For BASES 




H 


M 


Sec 


Overhead 













.12 


Grid Dptim. 


Step 





13 


52 


.62 


Integration 


Step 


1 


5 


42 


.20 


Go time for 


all 


1 


19 


34 


.94 



(2) Expected event generation time 
Expected time for 1000 events : 



> BASE1 

> BASE3 : 
BASE4 : 
BASE5 : 
RANDM : 
PLOTH : 
PLOTB : 
BSRSLT: 



0.40 Sec 



Making bases. rz. 



finished 
finished 
finished 
finished 
finished 
finished 
finished 
finished 



===> Directory : //bn 

3 (N) BASES 

1 (N) BASES 

2 (N) BASES 



data(real*8) 

data(integer*4) 

data(real*4) 
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A part of the standard output from spring 

Suppressing decay of the following particles in PYTHIA, 

K_S0 eta LambdaO Sigma+ SigmaO Sigma- 

Xi- Omega- D+ DO D_s+ Lambda_c+ 

tau- pi+ K+ K_L0 piO 



XiO 
mu- 



Loading 


bases . 


rz. . . 


> 


BASE1 


: finished 


> 


BASE3 


: finished 


> 


BASE4 


: finished 


> 


BASE5 


: finished 


> 


RANDM 


: finished 


> 


PLOTH 


: finished 


> 


PLOTB 


: finished 


> 


BSRSLT 


: finished 



** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 



Welcome to the Lund Monte Carlo! 



-><- 



ep 



e+e- 



PP 



PPP 
P P 
PPP 
P 
P 



Y Y 
Y Y 
Y 
Y 
Y 



T 
T 
T 
T 



H H 
H H 
HHHHH 
H H 
H H 



III 

I 

I 

I 
III 



A 
A A 
AAAAA 
A A 
A A 



This is PYTHIA version 6.136 
Last date of change: 30 Nov 1999 



Disclaimer: this program comes 
without any guarantees . Beware 
of errors and use common sense 
when interpreting results. 

Copyright T. Sjostrand (1999) 



An archive of program versions and documentation is found on the web: 
http : //www . thep . lu . se/~torb j orn/Pythia . html 



When you cite this program, currently the official reference is 
T. Sjostrand, Computer Physics Commun. 82 (1994) 74. 

** The supersymmetry extensions are described in 

** S. Mrenna, Computer Physics Commun. 101 (1997) 232 

** Also remember that the program, to a large extent, represents original 

** physics research. Other publications of special relevance to your 

** studies may therefore deserve separate mention. 
** 

** Main author: Torbjorn Sjostrand; Department of Theoretical Physics 2, 

** Lund University, Solvegatan 14A, S-223 62 Lund, Sweden; 

** phone: + 46 - 46 - 222 48 16; e-mail: torbjorn@thep.lu.se 

** SUSY author: Stephen Mrenna, Physics Department, UC Davis, 

** One Shields Avenue, Davis, CA 95616, USA; 

** phone: + 1 - 530 - 752 - 2661; e-mail: mrenna@physics.ucdavis.edu 

** 

** 



** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
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I****************** PYINIT: initialization of PYTHIA routines ***************** 



PYTHIA will be initialized for p+ on e+ user configuration 



P+ 
e+ 



px (GeV/c) 
. 000 
0.000 



py (GeV/c) 
. 000 
0.000 



pz (GeV/c) 
820.000 
-27.520 



E (GeV) 
820.001 
27.520 



corresponding to 300.444 GeV center-of -mass energy 



******** PYMAXI : summary of differential cross-section maximum search ******** 



I ISUB Subprocess name 
I 



I I 
I Maximum value I 
I I 



I 308 e+ uu~dd~ss~ -> e+ q m+ m- 
I 



5.1962D-02 I 
I 



######*###*###******** PYINIT: initialization completed ********************** 



=> START of SPRING at 1:21(13/ 4/ 0) 



Number of generated events = 100 

Event listing (summary) 



I particle/jet KS 


KF 


orig 


P- 


_x 


P- 


-y 


P- 


_z 


E 


m 


1 !p+! 


21 


2212 





0. 


,000 





,000 


820 


.000 


820. 


.001 


0. 


.938 


2 !e+! 


21 


-11 





0. 


.000 


0. 


,000 


-27 


.520 


27. 


.520 


0. 


.001 


3 !d! 


21 


1 


1 


-0. 


.343 


-0 


,058 


168 


.409 


168. 


.409 


0. 


.000 


4 !e+! 


21 


-11 


2 





,000 


0. 


.000 


-27 


.520 


27 


.520 


0. 


.000 


5 !d! 


21 


1 


3 


1. 


.167 


-0. 


.401 


138 


.481 


138. 


.487 


0. 


.000 


6 !e+! 


21 


-11 


4 





,000 


0. 


.000 


-27 


.497 


27. 


.497 


0. 


.000 


7 !d! 


21 


1 





-10. 


,149 


16. 


.228 


96 


.119 


98. 


.007 


0. 


.004 


8 !e+! 


21 


-11 





0. 


.009 


0. 


,020 


-6 


.958 


6. 


.958 


0. 


.001 


9 !mu-! 


21 


13 





19. 


.804 


-7. 


.523 


-12 


.433 


24. 


.564 


0. 


.106 


10 !mu+! 


21 


-13 





-8. 


,498 


-9. 


.126 


34 


.256 


36. 


.455 


0. 


.106 


11 e+ 


1 


-11 


8 


0. 


.009 


0. 


.020 


-6 


.958 


6. 


958 


0. 


.001 


12 mu- 


1 


13 


9 


19. 


,804 


-7. 


,523 


-12 


.433 


24. 


.564 


0. 


.106 


13 gamma 


1 


22 


4 


0. 


.000 


0. 


.000 


0. 


,000 


0. 


.000 


0. 


.000 


14 mu+ 


1 


-13 


10 


-7. 


,882 


-8. 


,464 


31 


.770 


33. 


.810 


0. 


.106 


15 gamma 


1 


22 


2 


0. 


,000 


0. 


,000 





,000 


0, 


.000 


0, 


.000 


16 (d) 


A 12 


1 


7 


-9. 


,379 


15. 


,009 


77 


.561 


79, 


.555 


0, 


.004 


17 (g) 


I 12 


21 


7 


-1. 


,386 





,556 


21 


.044 


21, 


.097 


0, 


.000 


18 (g) 


I 12 


21 


3 


-1. 


,510 


0. 


.343 


29 


,905 


29. 


.945 


0, 


.000 


19 (uu_l) 


V 11 


2203 


1 


0. 


,343 


0. 


.058 


651 


.591 


651. 


.591 


0, 


.771 



27 



20 


(string) 


11 


92 


16 


-11. 


.932 


15. 


,967 


780. 


.101 


782. 


,188 


53. 


,513 


21 


(rho-) 


11 


-213 


20 


-2. 


.614 


3. 


,803 


19. 


.032 


19. 


,599 


0. 


,763 


22 


eta 


1 


221 


20 


-5. 


,726 


8. 


,838 


49, 


.702 


50. 


.808 


0, 


.547 


23 


(K*+) 


11 


323 


20 


-0. 


,313 


0. 


,722 


5, 


.628 


5, 


.751 


0, 


.885 


24 


(K*-) 


11 


-323 


20 


-1. 


,601 


1, 


,458 


13. 


.870 


14, 


.064 


0, 


.862 


25 


(rho+) 


11 


213 


20 


-0. 


,642 





,548 


18 


.730 


18. 


.778 


1, 


.033 


26 


piO 


1 


111 


20 


-0. 


,205 


0. 


.308 


3. 


.644 


3, 


.665 


0, 


. 135 


27 


(rho-) 


11 


-213 


20 


-0. 


,597 


-0. 


.035 


42. 


.069 


42. 


.081 


0, 


.796 


28 


(rhoO) 


11 


113 


20 


-0. 


. 171 


0. 


.378 


46. 


.359 


46. 


.369 


0. 


.877 


29 


piO 


1 


111 


20 


0. 


.093 


-0. 


,222 


22. 


.348 


22. 


,349 


0. 


.135 


30 


K+ 


1 


321 


20 


-0. 


.371 





.069 


95. 


.962 


95. 


,964 


0. 


,494 


31 


SigmaO 


1 


3212 


20 


0. 


.358 


0. 


.006 


317. 


.588 


317. 


,590 


1. 


,193 


32 




1 


211 


20 


-0 . 


, 143 


. 


. 094 


145 . 


. 170 


145 . 


. 170 


. 


, 140 


33 


pi- 


1 


-211 


21 


-1. 


,237 


1. 


.578 


6. 


.910 


7. 


.196 


0. 


,140 


34 


piO 


1 


111 


21 


-1. 


.376 


2. 


.225 


12. 


.122 


12. 


,402 


0. 


.135 


35 


K+ 


1 


321 


23 


-0. 


.236 


0. 


.131 


2. 


.313 


2. 


.380 


0. 


,494 


36 


piO 


1 


111 


23 


-0. 


.078 





.591 


3. 


.315 


3. 


.370 


0. 


.135 


37 


(KbarO) 


11 


-311 


24 


-0. 


.822 





.679 


8. 


.481 


8. 


.562 


0. 


,498 


38 


pi- 


1 


-211 


24 


-0. 


.779 


0. 


.779 


5. 


.389 


5. 


.502 


0. 


,140 


39 


Pi* 


1 


211 


25 


-0. 


.837 


0. 


.227 


15. 


.119 


15 


,144 


0. 


,140 


40 


piO 


1 


111 


25 


0. 


,195 


0. 


.321 


3. 


.611 


3. 


.633 


0. 


.135 


41 


pi- 


1 


-211 


27 


-0. 


,260 


0. 


.270 


29 


.720 


29. 


,722 


0. 


,140 


42 


piO 


1 


111 


27 


-0. 


.337 


-0. 


.305 


12. 


.350 


12. 


,359 


0. 


,135 


43 


pi+ 


1 


211 


28 


-0. 


.149 


-0. 


,213 


8, 


213 


8, 


218 


0. 


.140 


44 


pi- 


1 


-211 


28 


-0 


.022 





,591 


38. 


.146 


38. 


,151 


0. 


.140 


45 


K_S0 


1 


310 


37 


-0. 


.822 


0. 


.679 


8. 


.481 


8, 


562 


0. 


,498 






sum: 


2.00 




0. 


.000 


0. 


.000 


792. 


.480 


847. 


.521 


300. 


.444 



===> Directory : //grp 

==========> END of SPRING at 1:21(13/ 4/ 0) 



1********* PYSTAT: Statistics on Number of Events and Cross-sections ********* 



I 
I 
I 

I 

I 

I N:o Type 
I 



Subprocess 



I 
I 
I 

-I- 
I 
I 
I 



Number of points 



Generated 



I 
I 
I 

1 

I 

Tried I 
I 



Sigma 
(pb) 



I I 

I All included subprocesses I 

I 308 e+ uu~dd~ss~ -> e+ q m+ m- I 

I I 



100 
100 



I I 

100 I 5.196D-02 I 

100 I 5.196D-02 I 

I I 



********* Fraction of events that fail fragmentation cuts = 0.00000 ********* 

************* << Total Cross-section >> ************* 

* * 

* ( 5.196186 +- 0.003313 )E-02 pb * 

* * 



28 



